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Abstract 

We present the first measurement of the electron angular distribution pa- 
rameter a2 in W ^ ev events produced in proton-antiproton collisions as 
a function of the W boson transverse momentum. Our analysis is based on 
data collected using the D0 detector during the 1994-1995 Fermilab Tevatron 
run. We compare our results with next-to-leading order perturbative QCD, 
which predicts an angular distribution of (1 it ai cos6* + 02 cos'^9*), where 6* 
is the polar angle of the electron in the Collins-Soper frame. In the presence 
of QCD corrections, the parameters qi and 02 become functions of p^, the 
W boson transverse momentum. This measurement provides a test of next- 
to-leading order QCD corrections which are a non-negligible contribution to 
the W boson mass measurement. 



I. INTRODUCTION 

After the discovery of the W boson [Q,^ at the CERN pp colhder, early studies of 
its properties verified its left-handed coupling to fermions and established it to be a spin 1 
particle [@,H. These were accomplished through the measurement of the angular distribution 
of the charged lepton from the W boson decay, a measurement ideally suited to pp colliders. 
The angular distribution was found to follow the well-known V — A form (l±cos6'*)^, where 
the polar angle 6* is the lepton direction in the rest frame of the W boson relative to the 
proton direction, and the sign is opposite that of the charge of the W boson or emitted 
lepton; this formulation assumes that only valence quarks participate in the interaction, 
otherwise the angular distribution is slightly modified. It is important to note that these 
measurements were performed on W bosons produced with almost no transverse momenta. 
This kinematic region is dominated by the production mechanism g -f g' — »• W. The center 
of mass energy used, y/s = 540 GeV, is not high enough for other processes to contribute 
substantially. 

At the higher energies of the Fermilab Tevatron {^/s = 1.8 TeV) and higher transverse 
momenta explored using the D0 detector ||^, other processes are kinematically allowed to 
occur. At low W boson transverse momentum, pj!^, the dominant higher order process 
involves initial state radiation of soft gluons. This process is calculated through the use 
of resummation techniques as discussed in Refs. 0-0. At higher values of pj^, where 



perturbation theory holds, other processes contribute [|T3|, such as: 

1. q + q' ^W + g 

2. q + g ^W + q' 

3. g + g^W + q + q' 

where only the first two contribute significantly at Tevatron energies |T^. These two pro- 
cesses change the form of the angular distribution of the emitted charged lepton to 

„ ; — OC (1 ± «! COS0* + Q;2C0S^ 6^*) (1.1) 

dpj;^ ay dcos 9* 



where the parameters ai and 02 depend on the W boson pt and rapidity, y |jT^. In Fig. |l|, 
the parameters ai and 0:2 are shown as functions of p^^. The angle 6* is measured in the 



Collins-Soper frame |]T5[; this is the rest frame of the W boson where the ^-axis bisects the 
angle formed by the proton momentum and the negative of the antiproton momentum with 
the X-axis along the direction of p^^. This frame is chosen since it reduces the ambiguity of 
the neutrino longitudinal momentum to a sign ambiguity on cos^*. 

In this paper, we present the first measurement of a2 as a function of pj^ |[T^, which 
serves as a probe of next-to-leading order quantum chromodynamics (NLO QCD), using the 
well- understood coupling between W bosons and fermions. This measurement probes the 
effect of QCD corrections on the spin structure of W boson production. 

At D0, the most precise W boson mass measurement is made by fitting the transverse 
mass distribution. However, since the transverse mass of the W boson is correlated with the 
decay angle of the lepton, the QCD effects discussed above introduce a systematic shift ~ 40 



MeV to the W boson mass measurement for events with p^^ < 15 GeV which must be taken 
into account. Presently, the Monte Carlo program used in the mass measurement models the 
angular distribution of the decay electron using the calculation of Mirkes |T^. During the 
next run of the Fermilab Tevatron collider (Run II), when the total error on the W boson 
mass will be reduced from the current 91 MeV for D0 
1 fb-i and to about 30 MeV for 10 fb~^ 



T7H22I to an estimated 50 MeV for 



23[| , a good understanding of this systematic shift 



is important. Therefore, a direct measurement of the electron angular decay distribution is 
important to minimize the systematic error. 

The paper is organized as follows: a brief description of the D0 detector is given in Sec. ^, 
with an emphasis on the components used in this analysis. Event selection is discussed in 
Sec. fll[ The analysis procedure is described in Sec. |^. Finally, conclusions are presented 
in Sec. |V|. 




100 125 150 175 
Pt^w [GeV] 

FIG. 1. The angular parameters ai (dashed) and a2 (solid) as functions of p^. These param- 
eters are evaluated integrated over the W boson rapidity, y. In the absence of QCD effects ai and 
02 equal 2.0 and 1.0, respectively. 



II. THE D0 DETECTOR 

A. Experimental Apparatus 

The D0 detector, described in more detail elsewhere 0, is composed of four major 
systems. The innermost of these is a non-magnetic tracker used in the reconstruction of 
charged particle tracks. The tracker is surrounded by central and forward uranium/liquid- 
argon sampling calorimeters. These calorimeters are used to identify electrons, photons, and 
hadronic jets, and to reconstruct their energies. The calorimeters are surrounded by a muon 
spectrometer which is composed of an iron-core toroidal magnet surrounded by drift tube 
chambers. The system is used in the identification of muons and the reconstruction of their 
momenta. To detect inelastic pp collisions for triggering, and to measure the luminosity, a 
set of scintillation counters is located in front of the forward calorimeters. For this analysis. 
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the relevant components are the tracking system and the calorimeters. We use a coordinate 
system where the polar angle 9 is measured relative to the proton beam direction z, and 
is the azimuthal angle. The pseudorapidity 77 is defined as — ln(tan|), and p is the 
perpendicular distance from the beam line. 

The structure of the calorimeter has been optimized to distinguish electrons and photons 
from hadrons, and to measure their energies. It is composed of three sections: the central 
calorimeter (CC), and two end calorimeters (EC). The r^-coverage for electrons used in 
this analysis is | r^ |< 1.1 in the CC and 1.5 <| r] |< 2.5 for the EC. The calorimeter 
is segmented longitudinally into two sections, the electromagnetic (EM) and the hadronic 
(HAD) calorimeters. The primary energy measurement needed in this analysis comes from 
the EM calorimeter, which is subdivided longitudinally into four layers (EM1-EM4). The 
hadronic calorimeter is subdivided longitudinally into four fine hadronic layers (FH1-FH4) 
and one course hadronic layer (CH). The first, second and fourth layers of the EM calorimeter 
are transversely divided into cells of size Arj x A0 = 0.1 x 0.1. The shower maximum occurs 
in the third layer, which is divided into finer units of 0.05 x 0.05 to improve the shower shape 
measurement. 



B. Trigger 

The D0 trigger is built of three levels, with each level applying increasingly more sophis- 
ticated selection criteria on an event. The lowest level trigger. Level 0, uses the scintillation 
counters in front of the forward calorimeters to signal the presence of an inelastic pp colli- 
sion. Data from the Level counters, the calorimeter and the muon chambers are sent to 
the Level 1 trigger, which allows the experiment to be triggered on total transverse energy, 
Eti missing transverse energy, $rp^ Ex of individual calorimeter towers, and/or the presence 
of a muon. These triggers operate in less than 3.5 /xs, the time between bunch crossings. A 
few calorimeter and muon triggers require additional time, which is provided by a Level 1.5 
trigger system. 

Candidate Level 1 (and 1.5) triggers initiate the Level 2 trigger system that consists of 
a farm of microprocessors. These microprocessors run pared-down versions of the off-line 
analysis code to select events based on physics requirements. Therefore, the experiment can 
be triggered on events that have characteristics of W bosons or other physics criteria. 

III. PARTICLE IDENTIFICATION AND DATA SELECTION 

This analysis relies on the D0 detector's ability to identify electrons and the undetected 
energy associated with neutrinos. The particle identification techniques employed are de- 
scribed in greater detail in Ref. p4|. The following sections provide a brief summary of the 



techniques used in this paper. 



A. Electron Identification 

Identification of electrons starts at tlie trigger level, wliere clusters of electromagnetic 
energy are selected. At Level 1, the trigger searches for EM calorimeter towers (A0 x At] = 
0.1 X 0.1) that exceed predefined thresholds. W boson triggers require that the energy 
deposited in a single EM calorimeter tower exceed 10 GeV. Those events that satisfy the 
Level 1 trigger are processed by the Level 2 filter. The trigger towers are combined with 
energy in the surrounding calorimeter cells within a window of A0 x A?] = 0.3 x 0.3. Events 
are selected at Level 2 if the transverse energy in this window exceeds 20 GeV. In addition 
to the Et requirement, the longitudinal and transverse shower shapes are required to match 
those expected for electromagnetic showers. The longitudinal shower shape is described by 
the fraction of the energy deposited in each of the four EM layers of the calorimeter. The 
transverse shower shape is characterized by the energy deposition patterns in the third EM 
layer. The difference between the energies in concentric regions covering 0.25 x 0.25 and 
0.15 X 0.15 in ?7 X must be consistent with that expected for an electron ^. 

In addition, at Level 2, the energy cluster isolation is required to satisfy fiso < 0.15, 
where fiso is defined as: 

_ Etotal(0-4)-£;EM(0.2) 

^-- £;em(0.2) ' ^'-'^ 

-Etotai(0.4) is the total energy, and i?EM(0.2) the electromagnetic energy, in cones oi R = 



(Arjy + (A0)^ = 0.4 and 0.2, respectively. This cut preferentially selects the isolated 
electrons expected from vector boson decay. 

Having selected events with isolated electromagnetic showers at the trigger level, a set of 
tighter cuts is imposed off-line to identify electrons, thereby reducing the background from 
QCD multijet events. The first step in identifying an electron is to build a cluster about 
the trigger tower using a nearest neighbor algorithm. As at the trigger level, the cluster is 
required to be isolated {fiso < 0.15). To increase the likelihood that the cluster is due to an 
electron and not a photon, a track from the central tracking system is required to point at its 
centroid. We extrapolate the track to the third EM layer in the calorimeter and calculate the 
distance between the extrapolated track and the cluster centroid in the azimuthal direction, 
pA0, and in the z-direction, Az. The cluster centroid position is extracted at the radius of 
the third EM layer of the calorimeter, p. The z position of the event vertex is defined by the 
line connecting the center of gravity calorimeter position of the electron and the center of 
gravity of its associated track in the central tracking system, extrapolated to the beamline. 
The electron Et is calculated using this vertex definition [E3]. The variable 



^.-i^Ui^f (3.2) 

where ap^ and cr^ are the respective track resolutions, quantifies the quality of the match. A 
cut of atrk < 5 is imposed on the data. Electromagnetic clusters that satisfy these criteria, 
referred to as "loose electrons," are then subjected to a 4-variable likelihood test previously 
used in the measurement of the top quark mass by the D0 collaboration p5[. The four 
variables are: 
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• A x^ comparison of the shower shape with the expected shape of an electromagnetic 
shower, computed using a 41-variable covariance matrix [^ of the energy depositions 
in the cells of the electromagnetic calorimeter and the event vertex. 

• The electromagnetic energy fraction, which is defined as the ratio of shower energy in 
the EM section of the calorimeter to the total EM energy plus the energy in the first 
hadronic section of the calorimeter. 

• A comparison of track position to cluster centroid position as defined in Eq. |3.2| . 

• The ionization, dE/dx, along the track, to reduce contamination from e'^e~ pairs due 
to photon conversions. This variable is effective in reducing the background from jets 
fragmenting into neutral pions which then decay into photon pairs. 

To a good approximation, these four variables are independent of each other for electron 
showers. Electrons that satisfy this additional cut are called "tight" electrons. 



B. Missing Energy 

The primary sources of missing energy in an event include the neutrinos that pass through 
the calorimeter undetected and the apparent energy imbalance due to calorimeter resolution. 
The energy imbalance is measured only in the transverse plane due to the unknown momenta 
of the particles escaping within the beam pipes. 

The missing transverse energy is calculated by taking the negative of the vector sum 
of the transverse energy in all of the calorimeter cells. This gives both the magnitude 
and direction of the J^rp, allowing the calculation of the transverse mass of the W boson 
candidates, M]f , given by 



Mf = ^2E^$t[1 - cos(0^ - 0*^)] (3.3) 

in which E^ is the transverse energy of the electron and (p^ and (p" are the azimuthal angles 
of the electron and neutrino, respectively. 



C. Event Selection 

The W boson data sample used in this analysis was collected during the 1994-1995 run of 
the Fermilab Tevatron collider. This data sample corresponds to an integrated luminosity of 
85.0 ±3.6 pb~^. Events are selected by requiring one tight electron in the central calorimeter 
(I r] |< 1.1) with Et > 25 GeV. The CC consists of 32 cj) modules. To avoid areas of 
reduced response between neighboring modules, the (p of an electron is required to be at least 
0.05 X 27r/32 radians away from the position of a module boundary. In addition, events are 
required to have ^^ > 25 GeV. If there is a second electron in the event (loose or tight) and 
the dielectron invariant mass Mge is close to the Z boson mass (75 GeV < Mge < 105 GeV), 
the event is rejected. 

To ensure a well-understood calorimeter response and to reduce luminosity-dependent 
effects, two additional requirements are imposed. The Main Ring component of the Tevatron 
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accelerator passes through the outer part of the hadronic calorimeter. Beam losses from the 
Main Ring can cause significant energy deposits in the calorimeter, resulting in false Ifirp. 
The largest losses occur when beam is injected into the Main Ring. Events occurring within 
a 400 ms window after injection are rejected, resulting in a 17% loss of data. Large beam 
losses can also occur when particles in the Main Ring pass through the D0 detector. Hence 
we reject events within a 1.6 /is window around these occurrences, resulting in a data loss 
of approximately 8%. After applying all of the described cuts, a total of 41173 W boson 
candidates is selected using electrons found in the central calorimeter. 

IV. EXPERIMENTAL METHOD 
A. Monte Carlo Simulation 

For this analysis, a Monte Carlo program with a parameterized detector simulation is 
used. This is the same Monte Carlo used in our previous results on the W boson mass 
measurement |]T^ and the inclusive cross sections of the W and Z bosons |^ , so it will only 
be briefly summarized here. 

In the Monte Carlo, the detector response is parameterized using the data from the 
experiment. This includes using Z bosons and their hadronic recoil to study the response 
and resolution. The response itself is then parameterized as a function of energy and angle. 

The kinematic variables for each W boson are generated using the RESBOS |1^ event 



generator with the theoretical model described in Refs. [|10| , p!3| , and the CTEQ4M parton 
distribution functions (pdf's) |2^. Finally, the angular distribution is generated according 
to the calculation of Mirkes f\A 



1. Hadronic Scale 

One of the parameters needed for the Monte Carlo program used in this study is the 
response of the calorimeter to the hadronic recoil, defined as the sum of all calorimeter 
cells excluding the cells belonging to the electron. The detector response and resolution for 
particles recoiling against a W boson should be the same as for particles recoiling against a 
Z boson. For Z — ;> ee events, we measure the transverse momentum of the Z boson from 
the e^e~ pair, p/*^, and from the recoil jet momentum, p-/*^'^, in the same manner as for 
W ^ ev events. By comparing prf:^ and P'^^'^., the recoil response is calibrated relative to 
the well- under stood electron response |]T^|. 



The recoil momentum is carried by many particles, mostly hadrons, with a wide momen- 
tum spectrum. Since the response of calorimeters to hadrons tends to be non-linear and the 
recoil particles are distributed over the entire calorimeter, including module boundaries with 
reduced response, we expect a momentum-dependent response function with values below 
unity. 

To measure the recoil response from our data, we use a sample of Z boson events with 
one electron in the CC and the second in the CC or the EC (CC/CC-I-EC). This allows 
the rapidity distribution of the Z bosons to approximate that of the W bosons where the 
neutrinos could be anywhere in the detector. Further, we require that both electrons satisfy 
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the tight electron criteria. This reduces the background for the topology where one electron 
is in the EC. We project the transverse momenta of the recoil and the Z boson onto the 
inner bisector of the electron directions (?7-axis), as shown in Fig. |^. By projecting the 
momenta onto an axis that is independent of any energy measurement, noise contributions 
to the momenta average to zero and do not bias the result. 

To determine the functional dependence of the recoil system with respect to the dielectron 
system, 'p,^^'^ ■ {—fi) is plotted as a function of p^?*^ ■ fj as shown in Fig. ^. For pf? > 10 GeV, 
the hadronic response is well described by a linear scale and offset: 



Pt ■ V = (^h Pt ■ V + PH 



(4.1) 



The parameters an and jSh are calculated using a least-squares fit to the data in the region 
p^T^ > 5 GeV, resulting in an = 0.972 ± 0.0095 and (3h = (-1.21 ± 0.14) GeV. For small 
values of p^, p^ < 10 GeV, the relation between the hadronic and electronic recoil is best 
described by a logarithmic function [|T^ ,p^: 



Pt'" ■V={1h ln(pr ' " v) + ^h) Pt'' ■ V 



(4.2) 



The parameters jh and 6h are derived using a least-squares fit to the data in the region 
pf < 10 GeV (see Fig. |), yielding jh = 0.099 ± 0.019 and Sh = 0.620 ± 0.047. In the 
intermediate region, 5 GeV < p^ < 10 GeV, the logarithmic and the linear fit match. 



P-^ 




p'rec A ; 



D ee A 



A 
11 



FIG. 2. Definition of the rj-^ coordinate system in a Z boson event, e/ denote the transverse 
momentum vectors of the two electrons. The r] axis is the bisector of the electrons in the transverse 
plane; the ^ axis is perpendicular to rj |19|. 



2. Tuning the Recoil Resolution Parameters 

In the Monte Carlo, we parameterize the calorimeter resolution, aj-^c, for the hard com- 
ponent of the recoil as 



Crec — SxecXi Pt 



(4.3) 
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50 60 70 80 

p^"»f^[GeV] 

FIG. 3. For Z — > ee events (points) the average value of p/®'^ • (—7?) is shown versus prp^ ■ fj. The 
Une shown is obtained from a hnear least squares fit to the data above p^ = 5 GeV as described 
in the text. The dotted lines represent the statistical uncertainties from the fit. 
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FIG. 4. For Z ^ ee events (points) the average value of p/^^ • (— ??) is shown versus prp^ ■ fj. 
Shown is the linear fit valid at p^ > 10 GeV and a logarithmic fit valid for p^ < 10 GeV. The 
dotted lines represent the statistical uncertainties from the linear fit. 



where s^ec is a tunable parameter, and p^*"^ is the recoil momentum of the hard component. 

The soft component of the recoil is modeled by the transverse momentum imbalance 

from minimum bias eventsQ This automatically models detector resolution and pile-up. 

To account for any possible difference between the underlying event in W boson events and 



^Minimum bias events are taken with a special trigger requiring only that a pp interaction has 
taken place. The kinematic properties of these events are independent of specific hard scattering 
processes and model detector resolution effects and pile- up which lead to finite ]^j^. 
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minimuin bias events, we multiply the minimum bias ^rp by a correction factor ttmb- We tune 
the two parameters Srec and amb by comparing the width of the ?7-balance, p-^^'^ -rj / Rrec+Pr^ 'V ^ 
measured from the CC/CC+EC Z boson data sample to Monte Carlo and adjusting the 
parameters in the Monte Carlo simultaneously until the widths agree. The width of the 
//-balance is a measure of the recoil momentum resolution. The recoil response, -Rrec is 
defined as 

R,^ = fell*!, (4.4) 

ml 

where qt is the generated transverse momentum of the Z boson. The contribution of the 
electron momentum resolution to the width of the ?7-balance is negligibly small. The con- 
tribution of the recoil momentum resolution grows with prf^ ■ fj while the contribution from 
the minimum bias ^j. is independent of prf^ ■ fj. This allows us to determine s^ec and amb 
simultaneously and without sensitivity to the electron resolution by comparing the width of 
the ?7-balance predicted by the Monte Carlo model with that observed in the data in bins 
of P'/'^ ■ fj. We perform a x^ fit comparing Monte Carlo and collider data. The values that 
minimize the x^ are found to be s^ec = 0.665 ± 0.062 GeV^/^ ^nd a^b = 1-095 ± 0.020. The 
non-linear hadronic scale in the region px <10 GeV leads to Srec = 0.50 ±0.06 GeV^/^, while 
Omb is unchanged. 



B. Extraction of the Lepton Angle 

Since only the transverse components of the neutrino momentum are measured, the 
transformation from the lab frame to the W boson rest frame (Collins-Soper frame) is not 
directly calculable. Therefore the polar angle of the electron from the W boson decay, 6*, 
is not directly measurable. In this analysis, 6* is inferred from the correlation between the 



transverse mass of the W boson and cos^* through the use of Bayes' Theorem p9 |. 

Experimentally, the only information we have about the W boson is that contained in 
the two kinematic variables M]f and p^^. But Mjf depends on the polar angle cos^*, the 
azimuthal angle 0* over which we have integrated, and pj!^. Therefore, the two experimen- 
tally measured variables M^f and p^^ give cos^*. An analytic expression exists for this 
relation (see Ref. 0]), so in principle the equation is solvable for cos 6**, but the experimen- 
tal values of both M^ and p^ include detector resolution effects that have to be unfolded to 
give the true cos 9* distribution. Even with perfect detector resolution, the equation would 
only be solvable if the W boson mass was known on an event by event basis. Therefore, 
we calculate the probability of measuring Mjf for a given value cos 6* in a given pj^ bin, 
p{M^\ cosd*,p^). This probability function is inverted to give the probability of measuring 
cos^* for a measured M^ , p{cos6*\M^ ,p^), using Bayes' Theorem: 

.fcosriM'^ p^) - P(^f|cosg%#Mcosr) 

p[cosu \Mt ,Pt)- jp(M]f|cos^*,p^Mcos^*)t/cos^* ^ ' 

where p(cos6'*) is the prior probability function, which we take as p(cos^*) = (1 + cos^ 6*), 
the charge-averaged expectation from V — A theory without QCD corrections. 
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To derive the probability function p{M^\ cos9*,pj!^), we use a Monte Carlo simulation of 
the D0 detector, which is described in Sec. [IV A . The correlation between Mjf and cos 6'* for 
p^ < 10 GeV is shown in Fig. ^. After determining j9(M7f| cos ^*,p^), it is inverted, yielding 
p{cos9*\M^ ,p^). The angular distribution is calculated by multiplying p{cos 9*\M^ , p^) 
with the measured transverse mass distribution. This is done in four p^ bins covering 0~10 
GeV, 10-20 GeV, 20-35 GeV, and 35-200 GeV. 

With the unfolded angular distributions now calculated, the value of a2 in each of the four 
p5^ bins can be determined. This is accomplished by generating a set of angular distribution 
templates for different values of a2- These templates are generated in a series of Monte 
Carlo experiments using the Monte Carlo program described in Sec. [IV A . 

The cos 6* templates are compared to the data through the use of a maximum likelihood 
method. Fig. |^ shows a series of angular distribution templates for different values of 02 
and p^ <10 GeV. 
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FIG. 5. Smeared W boson transverse mass versus true cos^* for p^ < 10 GeV from Monte 
Carlo. Acceptance cuts have been applied to events in this plot. This correlation plot is used to 
infer the cos 6* distribution from the measured M^ distribution. 



1. The Treatment of ai 



Since there is no magnetic field in the central charged particle tracking detector, it is not 
possible to identify the charge of the electron. Without charge identification, this analysis 
can only be performed by summing over the W boson charge and polarization. This implies 
that the linear term in cos 6* averages to zero in the limit of complete acceptance. However, 
after acceptance cuts have been applied, even the charge averaged angular distribution does 
depend on the linear term. The reason is that events generated with a non-zero ai correspond 
to slightly more central electrons after they are boosted into the lab frame compared to events 
generated with ai set to zero. After acceptance cuts have been applied, fewer events are 
lost at large cos^*. However, since this is only a second order effect, this measurement is 
not sensitive to cti. For this analysis, we calculate ai []T1| based on the measured p^^of 
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FIG. 6. Templates of the angular distribution for various Q2 values for p!p < 10 GeV. These 
templates are obtained from Monte Carlo after acceptance cuts have been applied which results in 
the drop-off at small angles. Each template is normalized to unity. 



each event. Possible variations of ai are treated as a source of systematic uncertainty (see 



Sec. IV E) 



C. Backgrounds 

To extract the electron angular distribution from the transverse mass distribution, the 
size of the backgrounds has to be estimated. The backgrounds are estimated as functions 
of the W boson transverse momentum and transverse mass, these being the two variables 
used to extract the angular distribution. The following sections describe how the four dom- 
inant backgrounds are calculated, and how they depend on transverse mass and transverse 
momentum. 



1. QCD 

A large potential source of background is due to QCD dijet events, where one jet is 
misidentified as an electron and the energy in the event is mismeasured resulting in large 
]^j.. This background is estimated using QCD multijet events from our data following the 
procedure described in detail in Ref. [|^. Briefly, the fraction of QCD background events 
in the W boson sample is given by 



fW _ ^1 
J QCD — AT 




(4.6) 



with the following variables: Ni and A^^ are the number of events in the W sample satisfying 
loose and tight electron criteria, respectively. The tight electron efficiency, e^, is the fraction 
of loose electrons passing tight cuts as found in a sample of Z boson events, where one 
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electron is required to pass tight electron identification cuts and the other serves as an 
unbiased probe for determining relative efficiencies. The jet efficiency, e^, is the fraction of 
loose "fake" electrons that pass tight electron cuts in a sample of multijet events. This sample 
is required to have low -^t(< 15 GeV) to minimize the number of W bosons in the sample. 
From this analysis, the overall QCD background fraction is found to be fqcD = (0.77±0.6)% 
with a transverse mass cut of 50 < M]f < 90 GeV imposed, this being the range used in 
the Bayesian analysis. For fg'cD ^s a function of p^, see Table |I[ 

2. Z^ee 

Another source of background is Z boson events in which one electron is lost in a region 
of the detector that is uninstrumented or one that has a lower electron finding efficiency such 
as that between the CC and the EC. This results in a momentum imbalance, with the event 
now being indistinguishable from a W boson event. This background can only be estimated 
using Monte Carlo Z boson events. The number of such Z boson events present in the W 
boson sample is calculated by applying the W boson selection cuts to herwig pi 



ee 



events that are processed through a geant [32| based simulation of the D0 detector and 



then overlaid with events from random pp crossings. This is done to simulate the underlying 
event, so that the effect of the luminosity can be included. The overall background fraction 
is found to be /^ = (0.50 ± 0.06)% averaged over all pj^. For the background fraction in 
each pj^ bin, see Table |I[ 

3. it Production 

The top quark background is not expected to contribute significantly, except in the 
highest p^ bin. The background from these events comes from t quarks decaying to W 
bosons. If one W boson decays electronically while the other decays into two hadronic 
jets, the event can mimic a high pt W boson event. This background, like the Z boson 
background, is calculated from Monte Carlo using herwig it events. The overall background 
fraction is f'^ = (0.087 ±0.027)%. For the background fraction in each p^^ bin, see Table |. 

4. W ^TV 

W ^>- TV events in which the r decays into an electron and two neutrinos are indistin- 
guishable from W ^ eu events. This background is estimated from Monte Carlo simulations 
using the W boson mass Monte Carlo described above. A fraction of the events is generated 
a.s W -^ Tu, decayed electronically, with acceptance and fiducial cuts applied to the decay 
electron in the same manner as in W ^ ev events. The acceptance for W —>■ tu —>■ evvu 
is reduced by the branching fraction B{t —>■ evv) = (17.81 ± 0.07)% ^3[. The kinematic 
acceptance is further reduced by the Et cut on the electron since the three-body decay of 
the r leads to a very soft electron Et spectrum compared to that from W ^ ei' events (see 
Fig. 1^). The fraction oi W ^ ru —>■ evvv events after these cuts are applied to the Monte 
Carlo is f^ = (2.03 ± 0.19)% over all pf. 
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For this analysis, the angular (cos6'*) templates are generated using the W boson mass 
Monte Carlo simulator with the branching ratio B{W -^ tv) = B{W ^>- eu), assuming 



lepton universality, and the above value for B{t - 
events (Fig. ^ is on average lower than that of W 
of the r. 



evv). The transverse mass oi W ^ ru 
> eu events, due to the three-body decay 




Electron E^ [GeV] 

FIG. 7. Electron Et spectrum for Monte Carlo W ^ tv ^ evvv events (dashed) and W 
events (solid histogram). Both spectra are normalized to unity for shape comparison. 



ev 




FIG. 8. Transverse mass distribution for W 
(dashed) from Monte Carlo. 



My^[GeV] 
ev events (solid) and W 



TV -^ evvv events 



5. Summary of Backgrounds 

As we have shown in the previous sections, and as can be clearly seen in Fig. ^, the 
background fractions in this measurement are small (a few per cent) over all M^ and p^ 
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ranges. The dominant backgrounds are due to QCD multijet events and Z boson decays, 
except in the highest p^ bin where the tt background is comparable in size. 




50 60 



90 



MtLGcV] 



70 80 90"" 50 60 70 80 

MtCGcV] MtI 

FIG. 9. Transverse mass spectrum for W ^ eu candidate events (solid histogram) and QCD 



(dashed), Z boson (dotted), and tt backgrounds (dashed-dotted) in four pj, bins. 



# [GeV] 


f^CD [%] 


f¥ [%] 


f¥ [%] 


0-10 


0.6 ±1.0 


0.16 ±0.02 


0.0028 ± 0.0009 


10-20 


1.0 ±1.0 


1.1 ±0.1 


0.025 ± 0.008 


20-35 


1.3 ±1.0 


1.4 ±0.2 


0.15 ±0.05 


35-200 


2.0 ±1.1 


1.7 ±0.2 


2.0 ±0.6 



TABLE I. Background fractions as a function of pj^ for events with a transverse mass cut of 
50 < M]f < 90 GeV imposed. 



D. The Measurement of a2 



To obtain the angular distribution for W boson events from data, the transverse mass 
distribution is inverted through the use of Bayes' Theorem as described in Sec. [IV B| . Since 
the probability distribution function used to invert the Mjf distribution is generated from 
Monte Carlo, we compare the background-subtracted Mjf distribution from data to that 
generated through our Monte Carlo to verify that it models the physics and detector correctly 
(see Fig. |l^). Based on a x^ test, the agreement between data and Monte Carlo is good; the 
X^-probabilities are 11.2%, 80.6%, 93.7% and 53.7% in order of increasing p^^ bins. Likewise, 
the experimental and Monte Carlo pj^ distributions can be compared, with the two showing 



agreement with a x^-probability of 7.4%, where only statistical errors are taken into account 
(see Fig. [n|). 

After extracting the angular distribution, the parameter 0:2 is computed using the method 
of maximum likelihood (see Fig. |l^). The angular distribution is compared to a series of 
Monte Carlo generated templates, each with a different value of 02 • The template that 
results in the maximum likelihood gives the value of 02 for each pj^ bin (Fig. |l^). The la 
uncertainties in 02 are approximately given by the points where the log-likelihood drops by 
0.5 units. To estimate the goodness of fit, the measured angular distributions are compared 
to these templates using a x^ test. The x^-probabilities that we obtain are 8.4%, 59.1%, 
87.7% and 11.6% in order of increasing p^ bins. 
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FIG. 10. Background subtracted transverse mass distributions (crosses) in four p^^ bins com- 
pared to Monte Carlo predictions (solid histograms). 



E. Systematic Errors 

Systematic errors on our measurement of 02 are due to uncertainties in the backgrounds 
and the parameters used to model the detector in the Monte Carlo. To estimate the errors 
due to the background uncertainties, the parameters from fits of the transverse mass distri- 
butions of the background are varied within their errors, and the analysis is repeated. For 
the errors due to detector modeling, the corresponding Monte Carlo parameters are varied 
within their errors and the analysis is repeated with new angular templates. For this analy- 
sis, we fixed ai to the values given by the next-to-leading order QCD prediction (see Fig. p. 
The error associated with this choice is estimated by changing ai to the value calculated in 
the absence of QCD effects (ai = 2.0). 
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FIG. 11. Background subtracted transverse momentum distribution (crosses) compared to 
Monte Carlo prediction (solid histogram). The error bars indicate statistical uncertainties only. 
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FIG. 12. Log-likelihood functions for four different pj^ bins. The arrows denote the values of 
maximum likelihood and the la errors. The vertical lines labeled V — A show ai = 1, the value 
for V — A theory without QCD corrections. 

The dominant systematic errors are due to uncertainties in the electromagnetic energy 
scale and the QCD background. All systematic errors are summarized in Table |I|. The 
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FIG. 13. Angular distributions for data compared to Monte Carlo templates for four different 
Pjf bins. Shown are the templates that fit best (solid) and the templates for q;2 = 1 (dashed) and 
02 = (dotted). 

systematic errors are combined in quadrature. The statistical uncertainties are, except for 
the first p}^ bin, larger by a factor of three than the systematic uncertainties. 



p^r [GeV] 


0-10 


10-20 


20-35 


35 - 200 


02, measured 


1.09 


0.84 


0.52 


0.13 


stat. errors 


±0.13 


±0.25 


±0.36 


±0.38 


02, predicted 


0.98 


0.89 


0.68 


0.24 


mean p}p 


5.3 


13.3 


25.7 


52.9 


QCD 


±0.04 


±0.05 


±0.09 


±0.07 


Z ^ee 


±0.01 


±0.02 


±0.02 


±0.04 


it 


±0.00 


±0.00 


±0.00 


±0.02 


EM scale 


±0.06 


±0.05 


±0.03 


±0.04 


hadronic scale 


±0.03 


±0.01 


±0.04 


±0.04 


hadronic resol. 


±0.02 


±0.02 


±0.05 


±0.06 


fixed oi 


±0.01 


±0.05 


±0.03 


±0.03 


combined syst. 


±0.08 


±0.09 


±0.12 


±0.12 



TABLE II. Central values for 02 with statistical and systematic errors. 
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F. Results and Sensitivity 



To estimate the sensitivity of this experiment, the x^ of the a^ distribution is calculated 
with respect to the prediction of the V — A theory modified by next-to-leading order QCD 
and that of the V — A theory in the absence of QCD corrections. The x^ with respect to the 
QCD prediction is 0.8 for 4 degrees of freedom, which corresponds to a probability of 94%. 
The x^ with respect to pure V^ — A is 7.0 for 4 degrees of freedom, which corresponds to 14% 
probability. To make a more quantitative estimate of how much better V — A modified by 
next-to-leading order QCD agrees over pure V — A, we use the odds-ratio method |, which 
prefers the former over the latter theory by ~ 2.3 a. The results of our measurement along 
with the theoretical prediction are given in Fig. and Table ||. 
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FIG. 14. Measured a2 as a function of pj^ compared to the next-to-leading order QCD calcu- 
lation by Mirkes (curve) and calculation in the absence of QCD (horizontal line). The combined 
systematic and statistical errors are shown as vertical bars, while the statistical errors alone are 
marked by horizontal ticks. 



V. CONCLUSIONS 

Using data taken with the D0 detector during the 1994-1995 Fermilab Tevatron collider 
run, we have presented a measurement of the angular distribution of decay electrons from 
W boson events. A next-to-leading order QCD calculation is preferred by ^ 2.3 o" over a 
calculation where no QCD effects are included. 



The odds-ratio R is defined as R 



n,P.(a2(NLOQCD)) , ,, , , 

H^ — -, — -, — ^^^-, where the product 
lliPi(a2(noQCD)) 



is over pjf bins. 



Pi(Q2(NL0 QCD)) is the normalized probability at the predicted value for 02 for the i*^ p^ bin. 



Pi (02 (no QCD)) is the normalized probability at the predicted value for V 
QCD effects, i.e. at 02 = 1.0. This corresponds to a la separation for log(-R) 



Pt 
A theory without 

0.5. 
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